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DE SITTER’S THEORY OF GALILEAN SATELLITES AND THE RELATED 

QUASI-PERIODIC ORBITS 


HENK BROER, LEI ZHAO 


Abstract. In this article, we investigate the mathematical part of De Sitter’s theory on the Galilean satellites, 
and further extend this theory by showing the existence of some quasi-periodic librating orbits by applications 
of KAM theorems. After showing the existence of De Sitter’s family of linearly stable periodic orbits in the 
Jupiter-Io-Europa-Ganymede model by averaging and reduction techniques in the Hamiltonian framework, we 
further discuss the possible extension of this theory to include a fourth satellite Callisto, and establish the 
existence of a set of positive measure of quasi-periodic librating orbits in both models for almost all choices of 
masses. 
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1. Introduction 

Based on the methods of studying celestial motions of planets by pure Keplerian approximations with 
additional corrections (called inequalities), formalizations of a mathematical theory of the Jovian Galilean 
satellites has been facing great difficulty for a long time. Indeed, even though this five-body system seems 
to be very much like a rescaled model of the inner solar system, and the mean motion, i.e. the Keplerian 
elliptic orbital frequency of Callisto is not close to lower-order resonance with the three inner satellites, the 
mean motions of lo, Europa and Ganymede are close to a 1 : 2 : 4 resonance, which give rise to inequalities 
greatly effecting the real behaviors of the system. 

In terms of averaging theory, we understand that the averaged system along this resonance differs by 
the corresponding (first-order) secular system, which is the averaged system of the perturbation over the 
three fast Keplerian motions, by some critical terms {i.e. resonant terms), corresponding to the 1 : 2 mean 
motion resonances of lo-Europa and Europa-Ganymede. The order of these resonance is so low such that 
these critical terms, in contrast to the critical term of high-order resonances, e.g. the 2 : 5 resonance roughly 
satisfied by the mean motions of Jupiter and Saturn, dominates, instead of being dominated by, the secular 
system. Moreover, not being first integrals of a first order approximating system (obtaining by adding the 
corresponding critical terms to the Keplerian part), the variations of the semi major axes, and therefore the 
contribution of the Keplerian part to the dynamics of the semi-fast variables near the resonances could also 
be significant. 

It is therefore important to analyze the dynamics of the first order approximating system obtained by 
adding up the corresponding critical terms, and to investigate whether it can be qualitatively continued 
when higher order effects are taken into account. These, together with a sufficient approximation of the 
normalizing transformation involved in the averaging procedure, provide a good theory for their celestial 
motions. This summarizes the strategy of W. de Sitter’s theory lldS251 . lldS31l of the Jovian Galilean 
satellites: lo, Europa, Ganymede and Callisto. 

As a starting point, De Sitter neglected Callisto and only considered a planar 4-body problem modeling 
the Jupiter-Io-Europa-Ganymede system. In modem terms, we may interpret this as follows: After being 
averaged over the 1:2: 4-resonance, truncated at the first orders of the eccentricities and of the masses, 
and further reduced by the rotational^ (9(2)-symmetry of this system, he established the existence of several 
families of periodic solutions with non-circular but almost-circular Keplerian ellipses, parametrized by one 
of the eccentricities of the satellites. Along these orbits, the semi major axes and the eccentricities of the 
satellites do not change. In a fixed reference frame, the pericentres (or “perijoves” when the massive body is 
called “Jupiter”) of the satellites process uniformly. De Sitter found this to be an important advantage of his 
theory ( lldS25l p. 8]). Indeed, in comparison with a theory built only on pure Keplerian motions, his first- 
order approximating system successfully removed the dynamical degeneracy of the uncoupled Keplerian 
systems that all their bounded orbits are closed, and very much dominates the local dynamics of the higher- 
order approximating systems and the full problem. Moreover, by a continuation method of Poincare, all 
these families of periodic orbits admit continuations to higher-order approximating systems and to the full 
system. Only one of these families of periodic solutions is found to be linear stable, to which, with properly 
assigned parameters, models and interprets the real motions of lo, Europa, and Ganymede in our solar 
system. 

Theorem 1.1. (De Sitter, ||dS09l ) In the planetary A-body problem with one mass sufficiently dominates 
the others, after symplectically reduced by the translation and rotational symmetry, there exists a family of 
linearly stable periodic orbits with 4:2: I-mean motion resonance. 
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Ganymede 



Figure 1. Jupiter-Io-Europa-Ganymede system at a collinearity 

In this article, we shall hrst represent these periodic orbits of the planar 4-body problem with some 
modern explanations and clarifications of this work lldS09l of De Sitter. For comparison purpose, we 
shall also treat what we would call 1 + 3-body problem of one hxed center and three satellites in parallel 
and establish the same result. Following De Sitter, we shall, after being reduced by the S 6>(2)-symmetry 
of rotations, first find particular periodic solutions of an approximating system, verify that they can be 
continued to periodic orbits of the full system, and identify the only linearly stable periodic solutions. 

In addition to this, for almost all masses, we shall establish the existence of librating quasi-periodic 
KAM orbits around the linearly stable family of periodic orbits. Since the real motions of lo, Europa, and 
Ganymede only roughly satisfy the 1:2: 4-resonance, these quasi-periodic orbits could serve to complete 
De Sitter’s theory of these satellites. 

Theorem 1.2. In both 1 -h 3-body problem or 4-body problem symplectically reduced by the translations 
and rotations, for almost all masses among which one sufficiently dominates the others, there exists a set 
of positive measure of quasi-periodic orbits librating around the family of linearly stable periodic orbits as 
indicated in Thm U.U 

In the analysis above, we have completely ignored Callisto, which is more distant from the Jupiter 
compared to the inner three satellites, but is nevertheless important to be included in the theory, since it 
might cause important secular inequalities of the inner three satellites. To have a more complete theory 
of the Galilean satellites including Callisto, we proceed by adding an outer fourth satellite to the system. 
Near circular motions of Callisto, De Sitter’s theory can be directly extended to such a model, in which 
we establish the existence of several families of quasi-periodic orbits, with some of them correspond to De 
Sitter’s periodic orbits, some of them correspond to librating quasi-periodic solutions established above. 

Theorem 1.3. In both 1 -h 4-body problem or 5-body problem symplectically reduced by the translations 
and rotations, for almost all masses among which one sufficiently dominates the others, there exist 

• a set of quasi-periodic orbits lying on a family of normally elliptic invariant 2-tori, along which the 
three inner mean motions satisfy the 4:2:1 resonance, and the outermost orbit is almost circular 
with frequency incommensurable with the frequencies of the inner three; and 

• a set of positive measure of quasi-periodic orbits, in which the three inner mean motions are close 
to 4 : 2 : 1 resonance, and the outermost orbit is almost circular with frequency incommensurable 
with the frequencies of the inner three. 
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We organize this article as follows: In Sectionj^we formulate our models and present them in canonical 
Joviancentric coordinates. A procedure of calculating the critical part of the perturbing function and the 
associated averaging procedure are presented in Section]^ Section [^presents De Sitter’s analysis on his 
family of linearly stable periodic orbits. We then present a version of KAM theorem in Section]^ which 
will be applied in Section]^ to establish librating quasi-periodic orbits around these periodic orbits, and in 
Section[7]to extend these results to the 1 - 1-4 / 5-body problems by establishing the existence of families of 
quasi-periodic orbits. 

We aim to present the main ideas with many more explanations in another expository work IBZI . which 
we hope could serve a broader community. 


2. Formulations: The 1 4- 3-Body Problem and the 4-Body Problem 

2.1. The 1 -H 3-Body Problem. To have a simplest model of our problem, we shall start by consider¬ 
ing the planar 1 -i- 3-body problem with a fixed center with mass mo and three small bodies with masses 
mi,m 2 ,m 3 << mo. In view of the fact that the mutual interactions of the bodies are smaller than the 
attracting forces from the center body, we may decompose the Hamiltonian F of the system as 


with 


F Kep — 


^ ~ ^Kep ^perti 


llpilP Ilp2ll^ llpsiP mm mmi mm 


nil m nij, roi ro 2 
m\m2 m2m2 m^rrii 


^03 


F — — 
^ pen ~ 


ri2 r23 

in whicfj^we have denoted the positions and momenta as 

{q\,q2,q3,p\,pi,pi) e T*({(^i,^2,?3) e 




X X ] 


^ ■ qi* qi,q2 + q 3 ,q 3 * q\]) 
and the mutual distances r,y = Wqi - qjW, 0 < i < j < 3, in which we have set qo = 0. 

2.2. The 4-Body Problem in Canonical Joviancentric Coordinates. In parallel, we shall also consider 
the planar full 4-body problem with masses mi, m 2 , m 3 << mo, positions qo, q\ , q2, q^ and momenta pQ,p\, P2,P3 
and write, following IILR95I . the Hamiltonian in the canonical Joviancentric coordinates 


(Pi), Pi, P2, P3, qi),qi,hy qi). 


defined by 

qo = qo 
qi=qi- qo 
§2 = q'2 - qo 
q 3 =q 3 - qo- 

The kinetic part thus becomes 


Po^ P0+ Pl+ P2 
Pi = Pi 
P2 = P2 
P3 = Pi- 


Pi 


T ^ To+ Ti 


with 


To^]:Y\\Pi 

^ i=l 


l"(- + —), 

rui mo 




0<i<j 


Pi-Pj . 
mo 


The potential part simply takes the form 


U ^Uo + Ui 


^Readers comparing this article with ldS09l should be careful that De Sitter took R = —Fp, 
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with 


'■oi' 


niim 


1=1 ‘ ^<i<j 

The complete Hamiltonian can thus be decomposed as the Keplerian part F^ep = Tq + Uq and the perturbing 
part Fpert = Ti + Ui- Compare with the 1 + 3-body formulation, the fact that Jupiter is not hxed causes the 
modihcation of the mass parameters in the Keplerian part, and the appearance of the indirect part Ti, which, 
though well-known to have no secular influence, does affect the study of lower order resonances (the reader 
is invited to consult the end of Section[^for this point). 

In this article, we shall treat both models altogether while commenting on their differences. 


2.3. Introduction of a Small Parameter; Order of the Mass Ratio. Let p be a small parameter repre¬ 
senting the order of the mass ratio between m,,; = 1,2,3 and mo. The momenta pi is of order yU if the 
velocities are considered as of order 1. We thus set 


mi = pmi,m 2 = ^im2,m-i = = ^ipi,p2 = PP2,Pi = m- 

Therefore, we have Fkcp ~ P and Fpen - 0{p^). 

We rescale the symplectic form by taking {pi,P 2 ,P 3 ,qi,q 2 ,q 3 ) instead of {pi,P 2 ,P 3 ,qi,q 2 ,q 3 ) as a set 
of Darboux coordinates. To preserve the same Hamiltonian vector field, we have to rescale the Hamiltonian 
function by a factor of yu"'. After the rescaling, we have (by abuse of notations) F^gp ~ 1 and Fpgg, - 0(p). 

2.4. Poincare’s Continuation Method. When p - 0, the system Fp^gp determines three uncoupled Keple¬ 
rian motions. This is a properly degenerate system. This proper degeneracy arise from the fact that in the 
Kepler problem all bounded orbits are closed, and is, very often, both a bless (the unperturbed dynamics 
is simple) and a difficulty (the perturbed dynamics of a degenerate systems cannot be well-controlled and 
higher order effects have to be considered to, hopefully, remove the degeneracy). 

In ||Poi92| . Poincare considered the problem of continuation of periodic orbits from proper-degenerate 
systems. 

Let (/, 6 ) - (/i, 01 , / 2 , 62 ) be a set of symplectic coordinates with (/i, 0i) e 7”*T"‘ and {I 2 , 02 ) e r*T"L 
F{1, 6 ) = ^ 0 ( 71 ) + p ■ Fiih = /«, e'i,l 2 , 62 ) + o(p) 

be a Hamiltonian function depending on a small parameter p. For small p, the dependence of pFi on Ii 
is much weaker than the dependence of Fq on these variables, therefore provided some non-degeneracy on 

Fq is satisfied, we could fix Ii in the expression of T^i. Along a periodic solution ^ of Fq, the frequencies 
dFo 

Hi - -are in complete resonance; we have denoted resonant angles (or critical arguments) by 9.. 

dll 

Theorem 2.1. ( Poincare. 0Poi921 n.46, p. 133]) The periodic solution ^ o/Fq can be continued to a periodic 
solution ofFforO < p « 1 provided that 

(P'Fo 

• The Hessian of Fq with respect to f is non-degenerate: det —r- + 0; 

dif 

• f corresponds to a non-degenerate critical point of Fi. 

In the original statement of Poincare, the second condition is weaker: It is enough to impose that f 
corresponds to a critical point of odd multiplicity of Fi. This “non-degenerate” version can eventually 
be phrased by the non-degeneracy of the associated Poincare map. The proof follows from the implicit 
function theorem. 
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2.5. First Order Approximating System and the Introduction of another Small Parameter: Order 
of the Eccentricities. Let us consider an invariant resonant 6-tori of Fkcp characterized by the 4:2:1- 
resonant condition on the mean motions. 

we average over this resonance in an Cst s-neighborhoocQ N (with a large 


To apply Thm 


2.1 


enough constant Cst) of these resonant tori, thus results in a certain function Fyes, containing only the critical 
terms, or resonant terms corresponds to this resonance, and truncate at the first order of the eccentricities. 
The function F is thus conjugated to 


FKep "t" Fyes "t" Fyeni, 

in which FKep is of order 1, Fyes is of order //e, while the remainder Fye^ is of the order 0 {fj.e^) -i- in 

which e is a small parameter characterizing the order of smallness of ei,e 2 ,ei. This averaging procedure 
will be made precise in Subsection |4.2| 


3. Calculation of the Critical Part of the Perturbing Function 


Let l\,l 2 ,h and g\,g 2 ,gi denote respectively the mean anomalies and the arguments of the pericentres 
of the particles qi,q 2 ,q 2 ,. We restricted ourselves to a deleted neighborhood of the circular motions to 
have these elements always well dehned. These elements, together with the semi major axes ai, 0 : 2.03 and 
ei, 62 , ^3 form a set of regular coordinates in such a deleted neighborhood in the phase space. 

To obtain the desired periodic solutions, we start by calculating the averaged perturbing function over 
the above-mentioned resonance, truncated at the first order of the eccentricities ei, 62 , ^3 that we suppose to 
be small. 

It is known (' 0Tis89l p.305]) that the function Fpeyi can be developed in Fourier series containing only 
cosines of the angles 

kill + k 2 lj + hgi + k4g j,{i, j) = (1,2), (2,3),(3,1),^ 1 ,^ 2 ,^ 3 ,^4 eZ 

with coefficients depending only on the semi major axes and the eccentricities. Among such terms we shall 
only be interested in those which are of first order in the eccentricities, and moreover, only those contain 
multiples of the critical arguments l\- 2/2 and 62 '-- h- 2 / 3 , since only these terms persist after being 
averaged over the resonance. 

As indicated in 0Tis89l p.307], to have terms at most of first order in the eccentricities, the corresponding 
{k\,k 2 ,k 2 ,k 4 ) must satisfy 

(1) |^i-^3l + fe-^4l< 1. 

Moreover, due to the invariance of the potential function under the rotations, we must have 

ko, + k^ — 0 . 

We see immediately that Eq. [1] cannot be satisfied for ^2 = 4^i + 0, i.e. the mutual interaction of the 
innermost and the outermost bodies q\,q 2 , could appear only at the second order of the eccentricities and 
are completely negligible for our present calculation of Fyes- For ^2 = 2^i, we can only have 

{h,k 2 ,k 2 ,k 4 ) = (l,-2,l,-l)or(l,-2,2, -2). 

It is thus sufficient to calculate the coefficients of each of the two corresponding terms (which are actually 

OT1OT2 OT2OT3 

of first order in the eccentricities) from the series expansions of-and-respectively. 

m f 23 

^We shall eventually only consider persistence of invariant objects under 0(e)-perturbations. The restriction is made to adapt to 


this. 
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mi m 2 

We now calculate, in the series expansions of-, the coefficients of the terms cos(li - 2/2 + ^1 - gi) 

r\2 

mi m 2 a\ 

and cos(li - 2/2 + 2gi - 2g2) while those in-are completely analogous. Let a - — be the semi 

ri2 02 

ri 

major axes ratio of the inner two bodies, p - — their ratio of radii, and vi, V 2 be the true anomalies of the 

r 2 

inner two bodies. Following IILR951 . we let 


( 2 ) — = - 


1 


'‘12 r2 VA + y’ 

in which 

A = 1 + - 2acos(Zi - ^2 + ^1 - g 2 ), 

V = 2a^cos(/i -^2+^1 -g2)- ^cos Z(qi,q2)j + P^ - a^- 

The function V is at least of first order in the eccentricities ei,e 2 , hence is much smaller than A in the 
deleted neighborhood we are now considering. We expand Eq. |^into 


(3) ± = 1 

ri2 r2 VA r2 VA^ 


+ 0(V^), 


The expansions of A 1,3 in terms of cos A:(li - I2 + .gi - ^2),^ = 0 , 1 , ■ ■ ■ is obtained by the usual 

Laplace coefficients bf\a), s — 1,1, defined as the coefficients of the Laurent series (c.f. BLR95II ) 


A ^ = (1 - az) ^(1 - az ^ ^ bf\a)z'‘. 

k=-oo 


A calculation shows that in first order of the eccentricities, r 2 only contains terms with argument I 2 , and 
that in first order of the eccentricities, V only contains terms with arguments 

h,h- gi + g2,2li - I 2 + g\ - g2, h, h + g\ - g2, h - 2 I 2 + gi - g2- 

We observe that to calculate these terms, it is enough to only calculate, in the expansions of A“'^^ and 
those terms with arguments 0, /i - ^2 + gi - g 2 , 2/i - 2/2 + 2gi - 2g2 and 3l\ - 3 I 2 + 3gi - 3g2 (c.f. IILR95I L 
The effective calculation is assisted by Maple 16. We find that 

^ “‘m 2 _ 21^ + 2gi - 2 g 2 ) -Be 2 cos(/i - 2/2 + gi - g 2 )} 


02 

{ - 


as 


(a 62 cos(/2 - 2/3 + 2g2 - 2g3) -Bei cosih - 2I3 + g2 - g3)} + c. 


- - a\ a2 

in which A, B are two functions of the semi major axes ratio a = — = —, 

a2 as 




3 

4‘ 

3 

4‘ 


1 


^ = 7«'^l/2(“) + 7«'^3/2(®) “ O' ^ 3 / 2 ( 0 ') - 7«^ 3 / 2 ( 0 ') 


and in which C is a certain function depending only on the masses and the semi major axes that will be 

ignored in the sequel. With a - -j=, we find, in tfie 1 + 3-body model, approximately 2A = 2.3810,2B = 

_ 

3.376^ The value of A is in consistence with De Sitter lldS091 for the 4-body model. The value of B differs 
from the 4-body model (for which De Sitter had noted the numerical value 2B - 0.964, while we found 

^Note that, in IdSO^ . De Sitter took (for small eccentricities) e,- = e,72, i = 1, 2,3 instead of in the expression of Fres- 
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2B - 0.8566 by Maple) by a difference of resulting from the indirect part of the perturbing 

function of the 4-body problem in the Joviancentric coordinates. 

The difference of the coefficient B in the two formalisms thus manifests the contribution of the indi¬ 
rect part of the perturbing function for the critical terms. To retain the simplicity of our analysis without 
diminishing its generality, in the sequel, we shall continue only with the assumptions A > 0, B > 0. 

4. Periodic Orbits of De Sitter 

From Thm |2.1[ we know that non-degenerate critical points give rise to periodic orbits of F for small 
fi. To look for such critical points, the particular form of Fres (that it contains only cosines) motivates the 
proposal of Poincare of the ansat 2 |^of collinearities of the three satellites when they simultaneously pass 
through their pericentres/apocentres (c.f. Figure [^. We therefore only consider those solutions of Fxep 
passing through those points with 

= { 0 ovTi ,0 ovn,Own, g2±n,g2,g2±T^), 

thus 32 families of them. These families are not entirely different. Indeed, since 6 i - l\ - 2 / 2 ,62 - h- 2/3 

does not vary along the resonance, leaving only the fast angle It, vary, the two families different only by tt 

in the fa-argument give rise to the same family of the periodic orbits. We therefore obtained 16 families of 
periodic orbits labelled by 

71 71 71 71 71 71 71 71 

(51,^2, ?71 := ^I -g2,m := g2 - gi) = (- ± - ± - ± - ± -)■ 

We name these families respectively by £)+ +_+,+, £)+_+,+ _ • • • etc., with subscript coincides with the ordered 
signs appearing in the ±’s of the right hand side. 

With this ansatz, a solution of F^ep + Fres is indeed periodic in some proper uniform-rotating frame, or 
in the reduced system by the S C>(2)-rotational symmetry, if and only if the frequencies of the precessions 
of the pericentres are equal. When such a periodic solution is found, its normal dynamics, as well as that 
of its continuation for small jj. can thus be determined by the corresponding evaluation of the Hessian of 

FKep Fres- 

4.1. Delaunay Coordinates and Reduction of the Rotational Symmetry. To effectively calculate the 
frequencies of the pericentres and to explicitly carry out the symplectic reduction of the rotation symme¬ 
try, we now introduce the canonical Delaunay coordinates (L,, f,, G,, g,,//,, h,),/ = 1,2,3, dehning as the 
following: 


Li = lUi ^fWi yjai 

circular angular momentum 

li 

mean anomaly 

Gi = Li ^1 - 

angular momentum 

gi 

argument of pericentre 


In which = m;,M, = mq in the 1 -H 3 -body problem, and = momi/(m{) + ^mi),Mi - niQ + ^fhi in the 
4 -body problem. 

To cope with our analysis near the 1 : 2 : 4 -resonance and to facilitate the reduction procedure, we mod¬ 
ify these Delaunay coordinates to a set of symplectic coordinates {D\,D2, DT,,di,d2, c/3, Zi,Z2,Z3,771,772,773) 

'‘in (dS03, De Sitter stated this as a necessity. This is not evident at all. 
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Figure 2. Planar Delaunay Angles 


in which 



II 

1 

D2 — 2Li + L2, 

62 - h- 2I3, 

7)3 = + 2Z>2 + ^3, 

^3 = h. 

N 

II 

0 

m = gi -82, 

Z2 - Gi + G2, 

^12 - 82- 83 , 

Z3 = Gi + G2 + G3, 

m = 83 - 


The coordinate Z 3 is just the total angular momentum; the rotational symmetry of the system translates into 
the independence of F of g^. The corresponding symplectic reduction procedure is thus achieved by fixing 
Z 3 and ignoring the variable gj. 


4.2. Elimination of the Fast Angle. We suppose 0 < oi < 02 < 03 ,0 < ei,e 2 , < e'^ < 1 such that 

ai(l+e^) ^ ^ 02(1 + ^ ^ 

02(1 - e ^) ’ 03(1 - e ^) 

so that the three elliptic orbits are bounded away from each other for all time. 

All these requirements determine an open subset P of the phase space. The coordinates 

(Di, D 2 , £> 3 ,^ 1 , 62 , ( 53 , Zi, Z 2 , Z 3 , 771 , t] 2 , r} 2 ) 

identifies to a subset of T® x R®. The variables 61,62 are critical (or semi-fast) arguments near the 
4:2: 1-resonance, with, near this resonance, the only fast angle < 53 . 

Let (Dj, D 2 , Dj) e be so chosen such that the set 

{Di =dO,D 2 = dO,D3 =d 0 ) 

consists in 4 : 2 : 1-resonant Keplerian motions in P. For Ci > 0, We define N to be the Ci//-neighborhood 
of the set {Di = DpD 2 = 0^,02 = Z)®) in P. We shall assume that Ci is chosen large enough to allow 
further application of KAM theorems. We let N to be the C 2 -neighborhood of {D\ - Dj, D 2 - D^, £>3 = D®) 
in P for some small enough C 2 so as to fixed a neighborhood for complex continuation of analytic functions. 
The small parameter fi is supposed to satisfy Ci/i < C 2 . 

Let Tq - X C® and {z & Tq : 3 z' & x R® s.t. \z - z'\ < s) be the s-neighborhood of 

T® X R® := R®/Z® X R® in Tc- Let Ta,.? be the s-neighborhood of a set A c T® x R® in T^. The complex 
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modulus of a transformation is the maximum of the complex moduli of its components. We use | ■ | to denote 
the modulus of either a function or a transformation. 

Being analytic functions on N, there exist i > 0, such that pKep and Fpen extend to analytic functions on 

T'n,s- 

For a function / : F® x R® ^ R, we define 



Proposition 4.1. There exists an 0(jS)-analytic symplectic transformation f : N ^ <PiN) such that 

f F — Fxep + F + Fygfji^ 

in which in N the analytic functions 

• Fres — {Fres}s, Is the truncation at the first order of the eccentricities of the function {Fpert)s,, ond 

• F re,n = 0{pe'^) + 

Proof We search for an auxiliary Hamiltonian H whose time-1 map gives the transformation f. We have 


(f) F - Fkcp + {Fpert + Xfj ■ FfCep) + P comp.V 

in which is seen as a derivation operator. Let Pres - {Ppen)ds be the average of Fpen over 6^, and 
Fpert = Fpert - {Fpen)ss be the zero-average part of Fpen- 

In the system pKep, since the frequencies VKep,\,VKep ,2 of 6\ and <52 is of order 0(p) compared to the 
frequency VKep,'}, of da, we do not have to solve the exact cohomological equation 

^Kep,\^Si^ ^Kep^l^S,^ ^Kep,3^ds^ ~ Fpert, 

instead, we just need H to solve the perturbed cohomological equation 

^Kep^i^ds^ — Ppert- 


We thus set 

1 ~ 

H - - I Fpert ddz 

^Kep,3 Jo 

as long as v/;ep ,3 + 0, which is indeed satisfied as a frequency of Keplerian elliptic motion. This amounts to 
proceed with a single frequency elimination for da. 

We have 

|H| < Cst p in 

We obtain by Cauchy inequality that in Tn,s-s„, \Xp\ < Cst \H\ < Cst p for some 0 < io < s/2. Shrinking 
from Tn,s- -So Tn>'.s-S0-s,, where N** is an open subset of N, so that fiTs-so-s,) Tn.s-So^ with s- So - 
Si > 0. The time-1 map f of thus satisfies \ f - Id\ < Cst p in Tn-,s-so-sc The function f*F is analytic 

in T^** s-so-s\ ■ 

The function F is thus analytically conjugated to 


(/> F — Fxep + Fres + Pcomp,\, 

and Wlemp 1 1 is of order O(p^): indeed, analogously as in | Fej02| , the complementary part 

PLmp.l ^ t)ft(^\ ■ pKep)dt -I- (jpiXp ■ FperJdt - VKep,2062^ “ VKep,3dd 


H 
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satisfies 


\Komp,l\ ^ Cst \X^\(\F^,,,\ + \Fpert I) + (yKep,2 + yKep,3)\H\ < Cst f/. 


To finish the proof, it is enough to expand F^es in powers of the eccentricities; F^es - Fres + 0(jj.e^), in 
which Fres only contains terms of hrst order in eccentricities. □ 


Remark 4.1. With the same method, the elimination procedure can be further continued to eliminate the 
dependence of ^3 in higher order perturbations as well (c.f. |Fej02|). 


The reduced system is thus determined by the Hamiltonian function 

FKep(D\,D2,Di) + Fres{Di, D2,02,61,62, Zi,Z2, 7 / 1 , J?!; Z3) + + o{f/) 

dehned on a subset of x R^, with Z 3 appearing as a parameter. 


4.3. Relative Frequencies of the Pericentres. We call the frequencies 

dFres 2 y/ 2 Acos( 6 i + 27 ji)e 2 m 2 + 2 .Bcos(bi + rji)eimi - y/2Acos(62 + 2772)617773 
dZi ylinae 162 

dFres 2 .Bcos(bi + 771)637771 - 2 V 2 Acos ((52 + 2772)637773 - 2 Bcos (62 + 772)637772 

V 2 = - = - ;— - 

dZ2 ■\lmoe2e3 

of 771 - gi - g 2 and ri 2 - g 2 - g 3 relative frequencies of the pericentres in Fxep + Fres- They appear as 
differences of the three quantities 

2V2Acos((5i + 2771)7772 2.Bcos(bi + 771)7771 - 2 '^Acos ((52 + 2772)7773 2 Bcos (62 + 772)7773 

^^moel ’ sJrnQe2 ’ sjme^ 

which are frequencies of g\,g 2 and g^ respectively. 

To obtain periodic solutions of Fkcp + Fres after reduction of the rotational symmetry, we have to ask 
= V 2 = 0. A simple analysis lldS091 pp. 10-12] with the signs of Vgi,Vg^,Vg^ shows that if we do 
not propose additional conditions on the masses, then this is only possible for families D- - +,+ (Case ( 6 ) of 
lldS09l l and D+ + +,+ (Case (16) of lldS09l ). Otherwise, with Q - '^Am 2 - 2 Bm\, we see that Q> 0, - - + 

(Case (2) of 108091), £>+.+,-,+ (Case (12) of |0SO9l), and for Q < 0, (Case (7) of |0SO9l), 

(Case (13) of lldS09l . see lldS091 p.lO]) are also possible for the nullity of vi and V 2 . In any of these cases, 
we obtain a one-parameter family of eccentricities ( 61 , 62 , 63 ) (parametrized by one of the eccentricities, e.g. 
by 62 ), such that the two relative frequencies vi, V 2 are zero. 


4.4. Continuation of Periodic Orbits. The existence of the corresponding families of periodic solutions 
of F/cep + Fres ^ftor reduction of the rotational symmetry thus follows from the above-mentioned families 
of eccentricities corresponding to the nullity of vi, V 2 ; The frequencies of 61 , 62 , 771,772 are all zero, leaving 
only the frequency of the angle ^3 non-zero. 

We now verify that Thm |2.1| is applicable in all these cases. The Hessian of F^ep with respect to Li, L 2 , T 3 
is a non-degenerate diagonal matrix. We therefore only need to verify that the Hessian of Fres with respect 
to (5i, 62 , Z\, Z 2 , 771 , 772 is non-degenerate. For this Hessian matrix, we summarize the results in the following 
table 
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Case 

Sign of Q 

Sign of the determinant of the Hessian matrix 


irrelevant 

positive 


irrelevant 

positive 


positive 

negative 


positive 

negative 


negative 

negative 


negative 

negative 


The additional perturbation is of the order Oip.e^) + Oip^). We thus obtain 

Proposition 4.2. The normalized mo, mi, m 2 , m 3 being fixed, after reduction of the rotational symmetry, for 
sufficiently small /i, e, each of families D+- + - of the periodic 

orbits of Ffcep + provided that their existence is permitted by the fixed value of m\ and fhj,, can be 
continued to families of periodic orbits ofF. 

Remark 4.2. In lldS09L the required non-degeneracy condition is not verified before the analysis of their 
stability. The non-degeneracy at the linearly stable periodic orbit is verified a posteriori by showing the 
non-degeneracy of its normal dynamics. 

4.5. Normal Dynamics/Linear Stability of the Periodic Orbits. We now determine the normal dynamics 
and the linear stability of these periodic orbits. For this purpose, an analysis for those periodic orbits in 
Fkcp + Fres suffices. For this, it is enough to fix D 3 and reduce the system by the S 6>(2)-symmetry of 
shifting (53, which reduce the periodic orbits to equilibria £’_ _ + + and respectively in the reduced 

system. 

We aim to calculate the eigenvalues of the linearization at the corresponding equilibra by using the 
Flessian of FKep + Fyes in coordinates {D\,D 2 , 6 \, 62 ,Z\,Z 2 , ? 7 i, 772 ) with small p. Since Fkcp and F^es do not 
appear in the same magnitude of p, it is enough to calculate the corresponding linearization matrix X from 
the Hessian Fl\ of Fkcp with respect to Di , D 2 and the Hessian FI 2 of Fres with respect to all other variables 
by 





X 

= J 

Diag{‘Hi,Fl 2 } 


in which 


'0 

0 

-1 

0 

0 

0 

0 

o' 



0 

0 

0 

-1 

0 

0 

0 

0 



1 

0 

0 

0 

0 

0 

0 

0 



0 

1 

0 

0 

0 

0 

0 

0 


J = 

0 

0 

0 

0 

0 

0 

-1 

0 



0 

0 

0 

0 

0 

0 

0 

-1 



0 

0 

0 

0 

1 

0 

0 

0 



,0 

0 

0 

0 

0 

1 

0 

0, 


We thus find that the linearization of the reduced flow at this point is a matrix of the form 
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'0 0 0 0 0 0 ' 

0 0 0 0 0 0 

yu-' yu-' 0 0 0 0 0 0 

yU-' yU-' 0 0 0 0 0 0 

0 0 yti^e 0 0 0 i?e 0 

0 0 0 0 0 0 

0 0 0 0 e-^ e-^ 0 0 

,0 0 0 0 e-^ e-^ 0 0 , 


in which all entries are only indicated by their orders in fi and e. 

We suppose that e is positive, non-zero and small enough. In this section, we aim to analyze whether 
the eigenvalues of this matrix is distinct and purely imaginary. For this purpose, it is enough to only argue 
with the small parameter fi and we may leave the parameter e aside for a moment. In Section]^ both of the 
orders of the frequencies in small parameters and e will be needed when we apply KAM theory to find 
librating quasi-periodic orbits. 

Due to the appearance of entries having negative powers of the small parameter yu, the calculation of the 
eigenvalues of this matrix is not convenient. This deficit is however found by De Sitter to be avoidable if 
we appropriately rescale the coordinates while preserving their symplecticity. 

We now treat 5,, 77 ,, i - 1,2 as periodic functions on R instead of functions on T, as otherwise we 
cannot continuously rescale them. Let cq - co(yu) be the rescaling factor to be determined. We take 
(Di,5i,D2,d2,Zi,77i,Z2, 772) as new coordinates with 


Di - 

D 2 - ■'^D2, 

Zi = 

Z 2 = 


= ^ 1 / 

^2 = ^ 2 / V^’ 
m = ml 
m = ml 


By choosing a posteriori cq - in these new coordinates, the corresponding linearization matrix 


takes the form (only the orders in yU 

are 

expressed) 






' 0 

0 

s/Ji 

0 

0 

0 

s/U 

0 ' 


0 

0 

0 

s/p 

0 

0 

0 

s/p 


s/p 

s/p 

0 

0 

0 

0 

0 

0 


s/p 

s/p 

0 

0 

0 

0 

0 

0 


0 

0 

s/Jl 

0 

0 

0 

s/U 

0 


0 

0 

0 

s/p 

0 

0 

0 

s/p 


0 

0 

0 

0 

s/p^ 

s[p^ 

0 

0 


. 0 

0 

0 

0 

s[p^ 

s[p^ 

0 

0 . 


with every entry containing a common factor This factor being ruled out, we arrive at a matrix having 
the form (again, only the orders in p are expressed) 
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'0 0 1 0 0 0 1 O' 

0 0 0 1 0 0 0 1 

1 1 0 0 0 0 0 0 

1 1 0 0 0 0 0 0 

0 0 1 0 0 0 1 0 ’ 

0 0 0 1 0 0 0 1 

OOOOyUjt/OO 

^0 0 0 0 jU 0 0^ 

we observe that by letting yU tend to zero, 4 eigenvalues of this matrix tends to 0, with 4 other eigenvalues 
tending to the eigenvalues of the upper-left 4 x 4-block, and thus their square powers satisfy the quadratic 
equation 

(4) - {Ans + A22 s’)x + (A 11 A 22 - A\2 )ss' - 0, 


in which 

. - 1 / 2 ^ pRep _i /2 ^ pRep - 1 / 2 ^ PRep 




dbl 


,S = —fi 


- 1/2 




Kep 


d 6 \ 


-,s = -{I 


- 1/2 




Kep 


dl\ 


dDydD 2 

The determination of the other 4 eigenvalues is achieved by a proper rescaling of the unknown variable 
A in the characteristic polynomial of this matrix. Indeed, in the corresponding determinant, if we replacing 
A by ■\/Jl A', we find that in the lower-right 4 x 4-block, the factor s/Jl appearing in the s/Ji A'’s can be ruled 
out by dividing out the factor in the last two rows and (then) the fifth and sixth columns. By letting ju 
tend to zero, we obtain a determinant, which (under the condition A 11 A 22 - Ajj > 0, ss' + 0, which could 
be verified easily) gives rise to a quadratic equation in (which is found to have the same form of Eq. |^; 


(5) 


• (Biicr -I- B22<j')x + (B 11 B 22 - B?2)o"0"" = 0’ 


in which 


D _ ..-m^^Pres o _ , - 3/2 9^ Pres 

t>i\ - B yPn - B - 


«9Zf 


5Zi5Z2 


,522 

dZl 


and 


- 1 / 2 ^ Pres - 1/2 


Of, 


(-■ 

\d 6 \df] 


;)/■ 


d 6 \ 


' g-' - ^-m ^^Pres ^ ^- 1/2 I d^Pres jd'^P, 


df]. 


\d 62 dfj 2 


}l- 


dh\ 


As such, we have obtained two quadratic equations. For small yu and small e, the corresponding periodic 
solutions are linear stable if both of these quadratic equations only have distincj^real, negative roots. 

The discriminant of Eq. I^is 


(AiiS H-A 22 S') - 4(AiiA22 -A\ 2 )ss' - (AiiS - A 22 S') +AA\ 2 Ss'. 


Since 


All < 0,A22 < 0 ,AiiA22 ~ Aj2 > 0, 

Eq. I^has two distinct real negative roots if and only if 


i > 0, i' > 0. 


As a result, by evaluating the values of s and s' at the corresponding points, we find that 
- are unstable. 

^Otherwise, property being not open, a lineai'ly stable periodic orbit is not necessarily continued to linearly stable periodic orbits 
for small pai'ameters. De Sitter seems to have overlooked this point, though it does not change the result. 
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Having verified that 811822 - 8^2 > 0 satisfied in all the other cases {i.e. for £_ _ +_+), 

the analysis of Eq. [^is similar. We thus find the following necessary and sufficient condition for Eq. [^to 
have two distinct real negative roots: 

cr > 0, cr' > 0. 

By evaluating the values of cr and cr' at the corresponding points, we thus hnd that £+ _ +_+, are 

unstable, and ^ ^ + is the only linearly stable family. In terms of periodic orbits, we have 

Proposition 4.3. (De Sitter, MdS09l ) is the only linearly stable family (parametrized by the small 

parameter p and the eccentricity e 2 ) of periodic orbits in the S 0(2)-reduced system of F^ep + Fres which 
can be continued to periodic orbits of F. 

Eor D- - + +, we hnd from the condition vi = V 2 = 0 that 

2 ■ 2 ^I^Am 2 e 2 sj 2 .Bm 2 e 2 

^1 — — i=z. -1—. S 3 = — 2=z. -1—• 

2 Y28mi + 2^1^Am^ 2 y28m\ + 2^1^Am^ 

The continuation of these linearly stable periodic orbits thus provide a mathematical explanation of the 
real evolution of the system Jupiter-Io-Europa-Ganymede, which was the base of De Sitter’s theory of the 
Galilean satellites lldS25l . 0dS31l . 

Question 4.1. Is the family _+_+ also linearly stable in the spatial problem? 


5. KAM THEOREM 


In order to establish the existence of librating quasi-periodic orbits, we shall make use of a general KAM 
theorem with parameters according to J. Eejoz | Eej04| , [ Fejl I) , | FejI0| . We note that a similar KAM 
theorem has also been established in I1BHT90II . 


5.1. Hypothetical Conjugacy Theorem. For p > 1 and q>0, consider the phase space R^xT'^xR^xR^ = 
{(/, 0, X, y)) endowed with the standard symplectic form dl A d 6 + dx A dy. All mappings are assumed to be 
analytic except when explicitly mentioned otherwise. 

Let 6 > 0,q' € {0,..., q}, q" - q - q', m e R^, and f e R®. Let be the (p + 2^)-dimensional closed 
ball with radius 6 centered at the origin in and - N^fi(d,q') be the space of Hamiltonians 

N e C'^iTP X Bf R) of the form 

c/ q 

N^c + {m,I) + Y^Pj(x]+y^^ + ^ /) + <A2(0),/®Z) + O^fZ), 

j=l j=q' + l 

with c e R, Ai e C"(T'’, R^ (8> R^), A 2 € C"(T^, R^ ® R^'*) and Z = (x,y). The isotropic torus T’ x {0) x {0) 
is an invariant tir-quasi-periodic torus of N, and its normal dynamics is elliptic, hyperbolic, or a mixture 
of both types, with Floquet exponents (3. The definitions of tensor operations can be found in e.g. I |F^ 

p.62]. 

Let y > 0 and f > p - 1, | ■ | be the ^^-norm on IP. Let HDy f = HDy Ap, q', q") be the set of vectors 
{'n 7 ,/ 3 ) satisfying the following homogeneous Diophantine conditions: 

\k-m + l'-fS'l^yUkf+ 1 )-^ 

for all k e ZP \ (0) and 1' 6 Z?' with \l[\ + ■ ■ ■ + |Z^,| < 2. We have denoted (J3i, ...,(3^,) by /?'. Let || ■ ||, 
be the s-analytic norm of an analytic function, i.e., the supremum norm of its analytic extension to the 
s-neighborhood of its (real) domain in the complexihed space x C.P jZP. 
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Theorem 5.1. Let {m°,fL') e HDy f and N° e j}"- For some d > 0 small enough, there exists s > 0 such 
that for every Hamiltonian N' e C"(T^ x such that 

||A^' - N‘’h < e, 

there exists a vector satisfying the following properties: 

• the map N' is of class C°° and is s-close to (w'’,yS'’) in the C°°-topology; 

• if (717, P) e HDy f, N' is symplectically smoothly conjugate to a Hamiltonian N e 
Moreover, e can be chosen of the form Csty^ (for some Cst > 0, k > 1) when y is small enough. 


Since analytic functions are C”, this theorem directly follow from the corresponding “Hypothetical 
Conjugacy Theorem” of | Fej04) . The fact that the aformentioned symplectic conjugation is actually analytic 
will appear in |Fej 111. 


5.2. An Iso-chronic KAM theorem. We now assume that the Hamiltonians N° - N° and N' - N( depend 
analytically on some parameter l e Recall that, for each i, N° is of the form 


'1 

N" = c” + «,/> + + Yi /> + <a, 2(0),/®z> + oalz). 


;=i 


J=9' + l 


Thm 5.1 can be applied to N° and N[ for each i. We will now add some non-degeneracy condition to the 
hypotheses of Thm 5.1 which ensures that 6 HDyf” actually occurs often in the set of parameters. 

Denote by 

//D" = [(m:,p:) e HDy^f : l e 


the set of “accessible” (y, rj-Diophantine unperturbed frequencies. The parameter is restricted to a smaller 
ball so as to avoid boundary problems. 


Corollary 5.1 (Iso-chronic KAM theorem). Assume the map 

is a dijfeomorphism onto its image. If e is small enough and if\\N( — N°\\d < b for each l, then for every 
(tir,y6) 6 HD^ there exists a unique l 6 B^^^ such that N[ is symplectically conjugate to some N € 
Moreover, there exists y > 0,f > p — such that the set 

[L e Bg : e //D“) 

has positive Lebesgue measure. 


Corollary 5.2 (Proper-degenerate Iso-chronic KAM theorem). Assume 

nO\ _ /^o,l /toJ „o,2 oO, 2 „0,h nO,h\ 

W.pJ -(tuy .mY .PY r" ,Pc ) 

such that for i j — dim{m°’\lfi''), — P + l ^tnc/ such that for any j — 1, 2 , • • • ,h, the frequencies involved 

in appear in the same magnitude of some small quantity ej. For small enough 1 >> ei » ■ ■ ■ » 

Cf, > 0 if the maps 

Vi : B'( R'A 11-> {m°'\ 
are dijfeomorphisms onto their images, then the map 

Vi : Bf ^ L ^ «,yS") 

is a dijfeomorphism onto their images. 

Consequently, we may choose s — Cste^y^ for large enough integer numbers N and k. 
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The proof goes in the same way as IIZhal4l Example-Condition 5.3], in which it is enough to replace 
the action variables by general parameters, and the Hessian matrix by the corresponding Jacobian matrix. 
Again, it is enough to observe that the determinant of the block diagonal matrix with blocks representing 
the Jacobian matrix of the mapping v,- dominates the rest terms in the expression of the determinant of the 
Jacobian matrix of v, as in the same spirit of OAKN061 Section 6.3.2]. 

6 . Invariant KAM Tori around the Linearly Stable Periodic Orbits of De Sitter 

Theoretically, a better mathematical theory for the Galilean satellites could be given by the quasi-periodic 
orbits on the possibly-existing invariant KAM tori around this linearly stable periodic orbit, since as noticed 
by De Sitter himself lldS09l . the 1:2:4 resonance is only satisfied roughly by the inner three Galilean 
satellites. In this section, we shall explore some extensions of De Sitter’s study by application of KAM 
theory. 

We have shown that for 0 < yu << 1 and 0 < e << 1, the elliptic equilibrium £_ _ + + of the SO(2) x 
S 6>(2)-reduced system Fxep + Fres is non-degenerate, with normal frequencies appearing at different orders 
in small quantities ji and e. 

We list the orders of these frequencies: 

(1) Frequency of the periodic orbit, which is the Keplerian frequency Vf^gp -^ of ^ 3 . It is non-zero and of 
order 1 . 

(2) Elliptic normal frequencies: It is direct to deduce from the dependence of the coefficients in Eqs. 
and|^ that forO<^(<<e<<l, two of them (denoted by Vn,\,Vn,i) are of order yfe, the other 
two of them (denoted by v„ 3, v„ i) are of order . 

With the hypothesis 0 < yu << e « 1, the two normal frequencies are small in their orders compared to 
the frequency of the periodic orbit: consequently v„,i, v „_2 is of smaller order compared to v„ 3 , Vn,A- 

To apply KAM theorems, we have to build higher order normal forms to control the smallness of the 
perturbation. We therefore consecutively eliminate the angle ^3 from the remainder F^g^ to have the higher- 
order remainder to be of the order G(yue^ -H yu^) for any prescribed large enough N by the same method as 
Prop |4.l| The system Ffcep + Fygs + Fggm is thus analytically conjugate to FKep + F^g^ + F^g^ in which 

• F’^g^ - Fggs = Oijie^ + yu^), 

• Fgg^ is independent of l^, 

• For N -2, Flg^ - Fres, and 

• F’^g,„ is of order 6 )(yue^) -H O(p^). 

We now analyze the dynamics of Fkcp + After being symplectically reduced by the rotational 
S (9(2)-symmetry, and the S (9(2)-symmetry of shifting I 3 , the reduced system is an (9(yU^ -Fyue^j-perturbation 
of the reduced system of FKep+Fres- Since the equilibrium point _+_+ of the latter is non-degenerate, this 
equilibrium point continues to exist for small enough yU and e, and give rise to an non-degenerate normally 
elliptic equilibrium + + of the reduced system of F^ep + which further give rise to a family of 
normally elliptic periodic orbits D^_ + + of F^gp + (with only the rotational S (9(2)-symmetry reduced). 
We thus reach the following 

Lemma 6.1. For any N and small enough yu and e, there exists a family of normally elliptic periodic orbits 
D^_ + + of Fxep + Fgg^ (reduced only by the rotational S (9(2)-symmetry) continuing the family D^_ + + of 

E Kep E rgs' 








18 


HENK BROER, LEI ZHAO 


We introduce another small parameter 6, which indicates the distance of a point to + + in the phase 
space of the fully reduced system. We develop the fully reduced system of F^ep + + + so that 

Fxep + = FKep + F^^^I + Oin^ed^) 

in which F^^^j is the corresponding quadratic part whose Hamiltonian flow is the linearized system at 
F^_ + It is thus of the order 

We are thus in a situation to perturb the integrable system Fkcp + F^^^ ^ by some perturbation of the order 
0 (p.^e^) + ( 9 (jU^) + 0 (j/e 6 ^). We impose that 


0 << ji « e « 1, and 0 < 5 << 1 small enough. 

The frequencies of those Lagrangian tori of F^ep + F ^^^, is dominated by the tangential-normal frequencies 
of E^_ + which are further dominated by the tangential-normal frequencies of E^_ + 

To apply Cor| 5 . 1 |, it is thus sufficient to verify the non-degeneracy condition for the tangential-normal 


frequencies Vper,v„^i,v„^2,yn,3,ynA of -,+,+ ■ Moreover, according to Cor 5.2 it is enough to verify the 
non-degeneracy conditions separately for frequencies in different scales. To avoid serious computational 
difficulties, we choose as parameters the rescaled masses mi,m2,m3, and the eccentricity 62- 


1 . Frequency of the periodic orbit; the required non-degeneracy condition is just 


dpKep 
dDi 


3t0. 


2 . Normal frequencies; This breaks down to the non-trivial dependence of the nontrivial coefficients of 
the monic quadratic equations with respect to parameters (the eccentricity 62 and the masses mi,m2,m3). 
We write Eqs|^and|^respectively as 


(6) 

and 


+ bix -I- Cl = 0 


( 7 ) 


+ b2X -I- C2 = 0 


respectively. 


Leituna 6.2. The Jacobians of (bi,ci) with respect to (mi, 62 ) and of (^ 2 ,^ 2 ) with respect to (m2, m3) are 
both non-degenerate almost everywhere. 


Proof. Assisted by Maple 16 , we And that 

evaluated at (mi = 1, m2 = m3 =0) equals to 

evaluated at (mi = l,m2 = m3 = 0) equals to -16 ■ 

The conclusion thus follows from analyticity and the fact that mi = /rmi, m2 = /im2, m3 = fiihi. 


d(bi,ci) 

d(mi,e2) 
d(b 2 , C 2 ) 

d(m2, m3) 


□ 


Therefore, for almost every fixed rescaled masses mi, m2, ™3 and fixed 62, there exists do > 0 , such that 
in a 5 o-neighborhood of the periodic orbits with fixed parameter 62 in the family Df_ ^ every Lagrangian 
tori of pKep + I has full torsion, in the sense that its frequency map (which is a small perturbation of the 
frequency map of the corresponding periodic orbit in the family D_ __++) is a local diffeomorphism. 


Theorem 6.1. Having fixed mo, there exists jUq > 0, such that for almost all masses mi, m 2 , m 3 satisfying 
max{mi,m2,m3) < jJ-omo, there exists a icf A,„, ,„2,m3 of positive measure consisting of invariant Lagrangian 
tori of F](gp + j close to the continued family D^_ ^ such that for any invariant Lagrangian torus in 
there exists a set of positive measure of masses pm'.^, pm'^, pm'^ with m'^,m 2 ,m'^ close respectively 
to mi, m 2 , m 3 , such that this invariant torus, with small deformation, persists under perturbation and give 
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rise to an invariant torus of the system F with mass parameter m', m'y These invariant tori of F form a 
set of positive measure in the direct product of the phase space with the space of masses {mi,m 2 , mf). 


Proof. Let (Ii,0i, 12 , 62 , 13 , 63 , 14 ,Of) be a set of action-angle coordinates of the fully reduced system of 
pRep + Pres I deleted neighborhood of such that for i - 1,2,3,4, the frequencies vP of 0, 

is close to the frequencies v„j respectively in the sense that v^. - v„_, -H 0(p). We localize ourselves in 
an ^-neighborhood of so that the action variables I\,l 2 ,l 3 ,h ~ d. In the system P/cep + Pj^^i 

only reduced by the rotational symmetry, such a neighborhood corresponds to a neighborhood of the family 
+ + of periodic orbits in which the invariant tori are obtained by fixing D 3 and /i, /a, /a, A. 

We see from Lem |6.2| that for almost all mi,m 2 ,m 3 , there exists an open set A in the phase space on 
which the map 

(D3, mi, m2, m3, ez) 1-^ (Vper, V«,l, V„, 2 , V„, 3 , V„, 4 ) 
is non-degenerate. Let p, e and ~6 be small enough such that the frequency map 


(D3, mi, m2, m3, 62) 1-^ 

is non-degenerate on an open set A°. Let T be any invariant torus in A". It is now enough to apply Thm |5.1| 
together with the help of Cor 5.2 to T for enough small parameters 0 < /t << e « 1, small enough 5 > 0, 
and large enough N. 

We now treat both the hypothetical frequency mapping tOi and the unperturbed frequency mapping u° 
as functions defined on the corresponding part of the direct product of the phase space with the space of 
masses (mi, m 2 , m 3 ). Since w, is C“ close to the unperturbed frequency mapping u°, their regular fibers 
have the same dimension. The last assertion thus follows from Cor l5.1l and the Fubini theorem. □ 


The following corollary follows from the above theorem and the Fubini theorem. 


Corollary 6.1. There exists po > 0, such that for almost all masses mi, m 2 , m 3 satisfying max{mi, m 2 , m 3 ) < 
Porno, there exists a set of positive measure of invariant Lagrangian tori of P in a small neighborhood of 
the continued family of linearly stable periodic orbits from the family + +. 

Remark 6.1. To verify the iso-chronic non-degeneracy condition near such a periodic orbit for a given set 
of masses, we shall have to calculate the corresponding Birkhoff invariants. Already with this relatively 
simple model, a direct calculation is hard to achieve even with the help of a computer. 


7. Extensions to the Planetary Five-body Problem 

We now add to the above-mentioned models a fifth far-away satellite “Callisto” with a small mass m^ 
and consider the corresponding 1 -F 4- and 5-body problems. We assume m^ ~ mi, m 2 , m 3 ~ p and write 
m4 = jum4. 

In the rescaled canonical Joviancentric coordinates {pi, p 2 , P 3 ,P 4 , qi, 62 , § 3 , § 4 ) (in which §4 designates 
the relative position of the fourth satellite with respect to Jupiter, p 4 is its linear momentum and pn - ppa), 
the Hamiltonian F of either of the systems (reduced from the translation symmetry in the five-body case) is 
analogously decomposed as 

P — P Kep P pert. 


Kep 


pert 


~ F + F XepA-f 

~ F pgpf + F pgpf/x. 


with 
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In which 


1 1/74 1P I^4M4 


^ Kep,4 


A^4 


ro4 


'' pert,4 


V P‘P* 
p/j — 


3 _ _ 

Z niiin/i 


To allow further analysis, we set 

Di^Lu 

Di — 2Ti + L 2 , 

7)3 = + 27,2 ^35 

A 4 = 7,4 
Zi =Gi, 

Z 2 -G\+ G 2 , 

Z' = Gi + G 2 + G 3 + G 4 , 

Z4 = L4 — G4, 

and, following Poincare, we further take 

^4 + m = V2(i4 


61 — h - 2 I 2 , 

62 — h - 2 / 3 , 

^3 = 

A4 - I4 + g4 - g3, 

r?i = ^1 -g2, 
m- § 2 - g3, 
m ^83- 

T]4 = -{g4 - 83)- 


■ G4)e 


-i(g4~S3) 


to have a set of regular symplectic coordinates in the neighborhood circular motions of the fourth body. 

We have thus obtained a set of Darboux coordinates 

(T>1, ^1, ^2> £>3> ^3> ^4> ^4> Zi, 7/1, Z2, 7/2, Z3,7/3, ^4, 774), 

which is regular up to circular orbit of the fourth body as long as the others are not. 

Recall that VKep ,3 is the Keplerian frequency of the third satellite. We denote by ^ the Keplerian 
frequency of the fourth satellite. 


With the hypothesis made in Subsection 4.2 we assume in addition that 03 < < 24,0 < 64 < < 1 such 


that 


773(1 + e^) 


< 1 


774(1 - e^) 

so that all the four elliptic orbits are bounded away from each other. The set (3 defined by these conditions 
is identified, by the above-mentioned Darboux coordinates, to a subset (3 of x x R^. 

Suppose that the set 

{Di =D?,D2 = D°,D3 =D°) 

consists in Keplerian motions in (3 for which the inner three frequencies satisfies the 4:2: 1-resonance. 
For Cl > 0, let f\A to be the transversal-Cantor subset of the Ci/ 7 -neighborhood of this set in <3 in which 
('^Kep, 3 , vkcpa) is (7> t)-Diophantine for some 0 < y << 1, t > 1. The fact that M. has positive measure for 
f > 1 and small enough y will be a consequence of the non-degeneracy condition of the Keplerian part with 
respect to D 3 and A 4 . This is part of the non-degeneracy conditions for application of KAM therem in the 
sequel, and is direct to verify by using the explicit formula of pKep- 

We shall assume that Ci is chosen large enough to allow further application of KAM theorems. We let 
(3 to be the C 2 -neighborhood of {Di - D^, D 2 — D^, D 3 = D®, A 4 = A®) in (3 for some small enough C 2 . 
The small parameter /7 is supposed to satisfy C1/7 < C 2 . 

Being analytic functions on Q, there exist s > 0, such that and Fpgri extend to analytic functions on 


Tas- 


For a function / : x R^ x R^ 


I, we define 


->277 


/ 7 / 537 //I 4 . 
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Proposition 7.1. For any integer N, there exists an 0(ju)-C°^ - Whitney symplectic transformation if/ : M ^ 
(p(M) such that 

if/ F — F + Ffgs + Ffgffi, 

in which in N the analytic functions 

• F fgs — {Fres}6j,A4- 

• Frem = 0{jie^ + 


Proof. The essential part of the proof is to consecutively eliminate the angles < 53 ,/I 4 with an elimination 
procedure analogous to the proof of Prop ITT] In Ad, under the Diophantine condition on the third and fourth 
frequencies, the existence of an C“-Whitney Hamiltonian solving the modified cohomological equation 
(dehned again by neglecting the (9(yu)-terms involving VKep.i and VKep.i) follows from e.g. | |Fg|T3] Prop 4], 
Finally, it is enough to truncate the resulting normal form at the (N - l)-th order of the eccentricities. □ 

Lemma 7.1. (Herman, see ||Fej04| Lemme 64]) The indirect part does not contribute to the 

’ -' niQ 

secular system in the sense that for i - 1 , 2 ,3 

j — Q 

Jv 'Wo 


The lemma is proven by observing that for i - 1,2,3,4, p, = while qi is a periodic function of 1, 
along Keplerian elliptic orbits, and application of the Fubini theorem. 

The Motion of Callisto is thus dominated by the corresponding secular system FsecA(^ 4 , ' 74 ), obtained 
from averaging the function 


yUmiOT4 fim2m4 

F pert,A — 

''14 ''24 


'•34 


over the fast angles ^3 and A 4 . We set 


F 


sec,A ~ 


I ^pert,4^^3^^4- 

Jt^ 


This is a function of order O(^e^) ( l|Tis89l p. 405]). In particular, it is of higher order compared to, and 
consequently dominated by F^es- We evaluate Acc ,4 at the corresponding circular orbits of the inner three 
to obtain a function Fsec. a- 


Claim 7.1. The point (0,0) is an elliptic equilibrium of Fsec.AifA, ^Ia)- 


This follows from the fact that the quadratic part of the analytic function Fsec.AifA, tIa) is even in (^ 4 , 774 ). 
More precisely, we deduce from 0Tis891 p. 405] that this quadratic part reads 

IJ.m4{mi]^B^'^\ai,a4) + m2^B^'^\a2,a4) + m3^B''Tfl3,fl4))(^4 + q^). 

in which ( irns89l p. 270, 271]) B^^\a, a') = 

We thus arrive at an approximating system: 

F Kep + Fees + F sec,A- 

We may thus identify the periodic solutions of the three satellites together with an circular orbit of the fourth 
satellite, which is used by De Sitter in ldS25ll as “intermediate orbits”. After symplectically reduced by the 
SO(2) X 5(9(2)-symmetry of shifting ( 3 , 44 , these solutions descends to an non-degenerate equilibrium of 
the reduced system, which can be continued to 

F Kep + F yes 
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for any N and small e,/i, with normal frequencies dominated by the normal frequencies of the corresponding 
equilibrium in the symplectically reduced system of F^ep + Fres + FsecA by this S 0(2) x S ( 9 ( 2 )-symmetry. 
Without reducing this S 0(2) x S ( 9 ( 2 )-symmetry, this equilibrium corresponds to normally elliptic invariant 
2 -tori of the system F^ep + F^es + FsecA "'bh an neighborhood by nearby librating invariant Lagrangian 
7 -tori. 

We consider the persistence of these invariant tori under small perturbation 0 (fie^ +fi^) for large enough 
N. We shall take the same parameters for the inner three satellites as in Section]^ The function FsecA(^4’ 74 ) 
contains the normalized mass OT4 as factor. To apply KAM theorem, we take 1714 and A4 as additional 


parameters. To see that the non-degeneracy condition is satished, in view of Cor 5.2 it is enough to notice 
further that this elliptic equilibrium is non-degenerate with normal frequency of order ju, and F^epA 
non-degenerate with respect to A4. 

We have thus proved the following theorems, by imposing 0 < fi « e « 1 and application of Thm [ 5 T| 
together with the help of Cor| 5 . 2 |in the same spirit as in the previous section: 


Theorem 7.1. When niQ is fixed, or almost all enough small masses mi, m2, m3, m4 and eccentricities e^, 64, 
the corresponding normally elliptic invariant 2-torus persists to exist for the system F with mass parameter 
m^,m 2 ,m'^,m 4 close to mi,m2,m3,m4, provided the small parameters p,e,6 are small enough. 


Theorem 7.2. When mo is fixed, for any fixed enough small masses mi, m2,m3,m4, there exists a set 
of positive measure consisting of invariant Lagrangian tori of F + F'f^^p such that for any 
invariant Lagrangian torus in Am,„, 2 ^ 3 ,m 4 . there exists a set of positive measure of masses m^,m' 2 ,m' 2 ,m'^ 
close respectively to mi,m2,m3,m4 such that this invariant torus, with small deformation, persists to ex¬ 
ist for the system F‘^ with mass parameter m^,m 2 ,m 2 ,m'^, provided the small parameters p,e,6 are small 
enough. These invariant KAM tori of F form a set of positive measure in the direct product of the phase 
space with the space of masses (mi, m2, m3, m4). 


We deduce the following corollary from the above theorem by Fubini theorem. 

Corollary 7.1. For almost all masses (mi, m2, m3, m4), there exists a set of positive measure of invariant 
Lagrangian tori of F in a small neighborhood of the continued normally elliptic invariant 2-tori. 
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